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ABSTRACT 
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Context. The Magellanic Cloud System (MCS) interacts via tidal and drag forces with the Milky Way galaxy. 
Aims. Using the Parkes Galactic All-Sky Survey (GASS) of atomic hydrogen we explore the role of drag on the evolution of the 
so-called Leading Arm (LA). 

Methods. We present a new image recognition algorithm that allows us to differentiate features within a 3-D data cube (longitude, 
latitude, radial velocity) and to parameterize individual coherent structures. We compiled an Hi object catalog of LA objects within 
an area of 70°x 85° (1.6 sr) of the LA region. This catalog* comprises information of location, column density, line width, shape and 
asymmetries of the individual LA objects above the 4-cr threshold of AT b 200 mK. 

Results. We present evidence of a fourth arm segment (LA4). For all LA objects we find an inverse correlation of velocities vqsr in 
Galactic Standard of Rest frame with Magellanic longitude. High-mass objects tend to have higher radial velocities than low-mass 
. ones. About 1/4 of all LA objects can be characterized as head-tail (HT) structures. 

Conclusions. Using image recognition with objective criteria, it is feasible to isolate most of LA emission from the diffuse Milky 
Way Hi gas. Some blended gas components (we estimate 5%) escape detection, but we find a total gas content of the LA that is about 
50% higher than previously assumed. These methods allow the deceleration of the LA clouds to be traced towards the Milky Way 
disk by drag forces. The derived velocity gradient strongly supports the assumption that the whole LA originates entirely in the Large 
Magellanic Cloud (LMC). LA4 is observed opposite to LAI, and we propose that both arms are related, spanning about 52kpc in 
space. HT structures trace drag forces even at tens of kpc altitudes above the Milky Way disk. 
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k> \ 1 . Introduction tion of the MS and LA. Not detecting of stars in the MS argues 

for drag forces, decelerating only the viscous gaseous compo- 

UO ■ Because of its proximity to the Milky Way, the MCS forms an nent? while the stdlar com p 0ne nt is unaffected and continues 

CTS ; ideal laboratory for studing the structure formation within the lo- its orbital motion The existence of th e LA, however still favors 

, cal universe in great detail. The more massive Milky Way galaxy strongly a tidal origin dNidever et a l. 2008). 
OO ' accrete the MCS and tidal forces and perhaps drag separate gas 

O from the stellar bodies. High-sensitivity optical studies disclose In thls P a P er we focus on the Hl 21 ~ cm emission distribution 

<N even the faintest stellar populations of evolved stars and allow of the LA as observed with the Galactic All-Sky Survey (GASS, 

■ to search for evidence of the stellar-gas-feedback processes iMcClure-Griffiths et alj2009t | Kalberla et al.| | 2010|). Assuming a 

K triggered by the intera ction with the Milky Way halo and/or lmear dlstance of about 50k P c from the Sun (|Cioni et al. || 2000|) , 

£ its gravitational field (llndu & Submmaniamll20ll . Both stel- single-dish radio telescopes like the Parkes 64-m dish offer a 

>< ' lar distributions are highly conc entrated, so only a fraction also linear resolution of about 200pc at Hi 21-cm wavelength. The 

H porjulate the Magellanic Bridge (Put man et alJMillrwin et al.1 hl S h s P ectral resolution of state-of-the-art radio spectrometers 

fl990h . So far the promin ent Hi structures denoted as the Leading not onl y allows determination of the bulk motion of individual 

Arm (LA. lPutman etZI IT998) and the Magellanic Stream have clouds but also examination of well resolved Hi line profiles. 

not detectable stella r counterpart (MS , Wannier & "wrixonll972t GASS for the first time makes full y sampled sensitive wide field 

iPutman et al.l[2003l:lBruns et al.ll2005h . ima S es of the Hl distribution available. 

Several investigation s of the MS make use of specialized sur- For our study we extracted a huge portion of the southern 

veys (Bri ms et ai]|2005l) or the Leiden- Argentine-Bonn all-sky sky from this data base to compile a comprehensive catalog of 

survey (Nidev er et al.ll2QQ8h . Eve n though a wealth o f informa- LA clouds and filaments (Sect.[2]). To achieve this aim we devel- 

tion has been compiled today (see lNidever et al.ll2008L and refer- oped an approach to differentiating between Milky Way gas and 

ences therein for a recent review), it is still a matter of debate to superposed (in space and frequency) MCS Hi emission. Based 

what extent tidal or drag forces determine the structure forma- on a newly developed extension of a standard image recogni- 
tion algorithm, we decompose the Hi data without the necessity 



* Table 2 is only available in electronic form at the CDS of Gaussian fits to the emission line profiles (Sect.[3). This al- 
via anonymous ftp to cdsarc.u-strasbg.fr (1 30.79.128.5) or via lows us to compile a comprehensive view of the LA down to a 
http://cdsweb.u-strasbg.fr/cgi-bonn/qcat 7J/A+A/ limit Tb = 200 mK. Using this catalog we investigate the ve- 
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locity structure of the LA and a newly identified leading arm 
four feature as a whole (Sect.|4]). The catalog comprises the ba- 
sic physical parameters of all identified LA objects; this gives us 
the opportunity to analyze the occurrence of head-tail structures 
and their orientation with respect to the MCS (Sect.0. We end 
the paper with a brief summary (Sect.©. 



to distinguish between different object classes according to their 
cataloged parameters, such as full area features (i.e. the Milky 
Way), "edge" objects (gas clouds touching the limits of the data 
cube either in position or velocity), small objects, and known 
galaxies. It is clear that the available source parameters must be 
complete enough to distinguish between different object classes. 



2. The data 

We use GASS data, observed with the L-band multi-feed 
receiver at the Parkes 64m telescope. Observations and 
initial data processing of t he survey were described by 



iMcClure-Griffiths et al] j2Q 09). We extracted data from the final 
data release (iKalberla et al.ll2010h that were corrected for stray 
radiation, instrumental effects, and radio-frequency interference 
(RFI). In calculating FITS cubes we used the web service 
as provided by htt p : // www.astro.uni - bonn.de I hi survey I but 
with two major extensions: 1) the positions were converted 
to the Mage ll anic S tream coordinate system as defined by 
iNidever et al.l (120081) . 2) The radial velocities were trans- 
for med into the galac tic- standard-of-rest system (GSR) (Eq. 2 
of lde Heii etaDl2QQ2h to quantify the motion of the MCS rel- 
ative to the Milky Way galaxy. We consider the GSR system 
as best suited to evaluating the dynamics of the coherent MCS 
structures relative to the Milky Way galaxy. This is most useful 
for visualization; however, it does not allow ambiguities to be 
overcome when LA clouds are blended with Milky Way gas. 

The GASS web server allows g enerating maps with differ- 
ent predefined angular resolutions (Kalberl a et al.l l201Qh . For 
our data analysis we chose two settings with effective beams 
of FWHM = 14.4arcmin (4.8arcmin per data cube pixel) and 
FWHM = 20.7arcmin, resulting in average rms uncertainties 
of 50 mK and 25 mK, respectively, for channel maps at a spec- 
tral resolution of The smoothed data were used to initialize our 
source finding algorithm (Sect. 13. IK all further processing was 
done on the high resolution data base. 

The whole field of interest covers 0° < /lmc ^ 70° and 
-42.5° < &lmc < 42.5°. This corresponds to 1.59 sr or 13% 
of the whole sky (Fig. [T] left). The radial velocity interval of in- 
terest covers -SOkms" 1 < vgsr ^ 250 km s" 1 . According to 
Ivan der Marel et al.l d2002), the LMC velocity is v G sr = 84 ± 
7kms _1 . Inspecting Fig. [2] (top panel) shows that the negative 
velocity limit cannot be determined very well because of the 
confusion with the Milky Way galaxy emission. The positive 
velocity limit safely includes all MCS features above our 4-cr 
threshold of T B = 200 mK (Fig. Q] right). 

3. Decomposing the LA emission 

Our aim was to compile an object catalog for the LA area, con- 
taining all data that are necessary for any interpretation, such 
as position, velocity, brightness temperature, center of emission, 
geometrical center, column density, and velocity gradient. For an 
accurate parameterization we also decided to determine first mo- 
ments in the position- velocity domain. This allows us to investi- 
gate the cloud's morphology and quantify their shape, kinemat- 
ics, and column density distribution asymetries in detail. The re- 
sulting catalog comprises emission features of all kind of struc- 
tures, regardless of their physical origin. 

A major task is to eliminate spurious features and to dis- 
tinguish between so urces belonging to LA and the Milky Way. 
INidever et al.l (120081) use a Gaussian decomposition for this pur- 
pose, which is a very elaborate enterprise. Intuitively, a param- 
eterization in first moments, as mentioned above, also enable us 



3.1. Image J and source identification 

The cloud catalog is automatically compiled by a newly de- 
veloped source finding and parametrization program within 
the ImageJ en yiroment. ImageJ is a scientific image process- 
ing software dRasbandl [1997-2009) that is most commonly 
used in biology and medicine but also partly in astron- 
omy. It is an open source software, and it can be found at 
http : //r sbweb . nih . gov/ij /index . html. 

The source-finding algorithm is designed to identify both 
spatially and kinetically coherent Hi objects above a given 
brightness temperature threshold. For each object some essen- 
tial physical parameters are derived, as presented above. 

3.1 .1 . Cleaning the catalog 

First we clean our database for insignificant features below the 
instrumental noise, for RFI, and residual systematical baseline 
defects. We use the smoothed data with a doubled kernel size of 
0.25°. All features below the threshold of T B = 200 mK can be 
excluded by this approach. Based on this threshold, we gener- 
ated a mask for the high-resolution data cube. 

The threshold mask isolates coherent structures and different 
individual objects. These structures need to be separated from 
other objects in phase space (spatially but also in velocity). This 
procedure attributes each pixel above the threshold to an asso- 
ciated object, permitting the calculation of the objects physical 
parameters. 

The GASS beam FWHM (14.4arcmin= 0.24°) corresponds 
to three pixels in the data cube. By this, any detected objects 
smaller than a single Parkes beam were excluded from further 
analysis. With this step we mitigate the imprint of RFI. We cross 
correlated our catalog with the HIPASS bight galaxy catalog 
(Koribalski et al. 2004). Within the area of interest, 361 galax- 
ies are located. Most of them are automatically excluded from 
further consideration because their velocities are beyond the ve- 
locity range considered here. Only three residual galaxies were 
needed to be accounted manually. 

Cold gas with an FWHM line width of up to three spec- 
tral channels (2.4 kms" 1 ) are considered to trace mainly the 
Milky Way emission because a sufficient high gas pressure is 
needed to warrant c old gas of ^ 250K in a two-phase medium 
dWolfire et al. 1995). HVCs containing extremely cold dense 
cores will still be cataloged because next to the narrow Hi emis- 
sion the enveloping warm gas will also be detected in adjacent 
spectral channels. RFI events which are however, usually nar- 
rowly confined in frequency, will be rejected. 

After such a selection, the LA catalog in its final version 
has 449 entries. In the following we show that these cataloged 
objects might safely be considered as belonging to the LA be- 
cause of their coherent velocity structure. Towards the field of 
interest, next to the LA, th e Wannier complexes are well known 
dWannier & Wrixonl 19721) . These HVC complexes (annotated in 
Fig. [2] and discussed in Sect. 4) form coherent structures but are 
separated from most of the LA structures by their lower GSR 
velocities. 
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3.1 .2. The Milky Way galaxy as an "edge" object 

Based on the catalog it is feasible to generate data cubes show- 
ing different representations of the LA objects. The spacial dis- 
tribution of all compiled objects is shown in Fig.QIleft panel). 
It includes objects that are boundary or edge objects of the data 
cube. These edge objects were excluded from further consid- 
eration (see FigQ] right panel). The comparison of both panels 
make it evident that the Milky Way disk emission is success- 
fully removed from the analyzed catalog. The Milky Way disk 
Hi emission above the threshold of 200 mK forms a single broad 
and coherent structure which touches all boundaries of the whole 
data cube. By definition the Milky Way galaxy Hi emission is 
classified as an "edge" object. Accordingly, when dropping the 
edge objects from the final catalog, the Milky Way emission is 
entirely removed. 

It is a big advantage of the applied algorithm to distinguish 
easily between the LA and the Milky Way disk without the need 
to fit other parameters. However, Hi objects physically belonging 
to the MCS but apparently "connected" with the Hi emission of 
the Milky Way galaxy will not be part of the final catalog. Here 
visual inspection might further improve the completeness of the 
catalog but could also introduce systematic biases. 



3.2. The cataloged parameters 

An extensive set of parameters is calculated for each cataloged 
object. The basic parameters are the position, the spatial and the 
kinematical centers, the intensity weighted center, the geomet- 
rical center, and the position of maximum brightness temper- 
ature. Further compiled parameters are the object extent in all 
three cube dimensions, the maximum brightness temperature, 
the maximum column density, the mean column density, the 
integrated brightness temperature, the phase space volume, the 
solid angle, the orientation angle, and the FWHM velocity width. 
The kinetic gas temperature is derived from the full- width-half 
maximum of the Hi line of the peak brightness temperature po- 
sition. Assuming a 50 kpc distance of the LA, the linear extent 
and the mass are calculated as well (see Tab.0. 

Figured displays numerous features that appear well sepa- 
rated from the diffuse Galactic emission. Only clouds are con- 
sidered that are unambiguously separated in phase space from 
the Milky Way. From a comparison of Fig. [T] (right) with Fig. [3] 
one might expect that some fraction of the total mass assigned to 
the leading arm might be missing. Correspondingly we find few 
objects at Magellanic longitude /ms -33° visible in Fig. [2] We 
use the area filling factor for the affected region in velocity and 
position (the projected phase space) as a measure to estimate the 
mass deficit, as less than 5%. 



3.3. Mass comparison 

The masses of LA features derived here are about 50% higher 
as recognized previously. In Table 1 masses of LA su bstructures 
are compiled and compared with iBriins et aD (120051) . The ob- 
vious differences are first be cause of different definition of the 
feature extents (in general in terims et al.l (120051) the features are 
systematically smaller). In particular LA 1.1 and LA 1.3 are ex- 
tended objects tha t are partly located outside the survey area of 
Brim s et al .(200!). 



Table 1. M ass comparison of the specific LA areas indicated in 
Fig. [T] with Briins et al . (2001). 



Region 


M ffl 
(1O 6 M ) 


M ffl (B05) 
(10 6 M ) 


Mm 
Mm (505) 


LA1.1 * 


2.67 


0.99 


2.69 


LAI. 2 


4.02 


3.14 


1.28 


LAI. 3 * 


6.38 


3.64 


1.75 


LAI * 


13.07 


8.26 


1.58 


LA2 


11.42 


9.09 


1.26 


LA3 


8.89 


7.44 


1.19 


LA4 


2.44 






Rest 


2.34 






LA total * 


38.15 


24.79 


1.54 



Notes. Masses from Briin s et al] (120051) are indicated by B05. All 
masses are calculated assuming a linear distance of 50 kpc. The asterisk 
(* ) mark areas that ar e obviously not congruent; i.e., the studied areas 
bv lBriins et al] {2005) are in general of smaller extent. 



4. The velocity structure of the LA as a whole 

We first study the velocity distribution of all LA objects to- 
gether to investigate all the environmental and orbital param- 
eters. To this aim we analyzed the radial velocity pattern of 
the cataloged objects in the GSR rest frame. Figure[3] gives an 
overview, while Fig. [2] shows details about the individual LA 
objects listed in the Appendix. To distinguish between LA fea- 
tures a nd Wannier complexes, w e included positions of Wannier 
clouds ( Wannier & Wrix onll 19721) in the upp er panel of Fig. [2] as 
given in the Wakker & van Woerden (1991) catalog. For /ms ~ 
60° WA and WB features cannot be distinguished easily from 
scattered LA features. LA2 and LA3, however, as is obvious 
from Fig. [3) stand out clearly. The complexes WC and WD are 
found at low GSR velocities, so separated well from LA clouds. 

The diameters of the circles in the lower panel are propor- 
tional to the third root of the cloud masses and specific LA sub- 
structures are colored. All together there is an obvious relation 
between /ms and vqsr- This finding implies that the cataloged LA 
as a whole has a coherent and continuous velocity structure. The 
clouds that are the most distant from the LMC tend to have the 
lowest velocities. Gas within the immediate vicinity of the LMC 
apparently had less time to get decelerated by the halo medium. 

At the location of the Galactic disk (around /ms = 30°, 
Fig.[2]bottom panel), there is a V-shaped underabundance of LA 
objects. We attribute this under- abundance to limitations of our 
approach to separate Milky Way Galaxy and LA at these posi- 
tions. HI clouds in this region have been attributed to Galactic 
emission. 

Connecting the location of the cataloged LA objects with 
their Hi mass (assuming D = 50 kpc) discloses a remarkable 
trend (in particular for LA 2 and LA 3): the higher the mass of 
the objects, the higher their radial velocities (Fig. [2) bottom panel 
and Fig.0middle panel). This argues for drag which decelerates 
the low-mass diffuse objects significantly while the high mass 
objects keep most of their momentum. In the immediate vicinity 
of the LMC (/ms < 25°, Fig. [5] top and bottom left panels), the 
higher mass objects do not behave any differently in velocity. 
Here, the low-mass objects are also observed at high vqsr radial 
velocities. Starting around /ms > 35°, the higher mass objects 
reveal that the highest radial velocites and lower mass objects 
are all at lower velocities. 
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70° 60° 50° 40° 30° 20° 10° 60° 50° 40° 30° 20° 10° 0° 

Magellanic longitude Magellanic longitude 



Fig. 1. Left: The map shows the column density distribution of the Hi 21 -cm line emission across the velocity interval -80 km s" 1 < 
vgsr ^ 250 km s" 1 . The complex Hi structure is produced by the superposition of the Milky Way galaxy emission, the LMC, and 
the LA. Right: The column density distribution of LA clouds across -80kms _1 < vgsr ^ 250 km s" 1 starting at a threshold of 
AT B = 200 mK (4 - <x level of the GASS data). All the objects shown are cataloged and parametrized. Coherent structures like the 
Milky Way disk and LMC emission are according to the applied criteria "border" objects and not part of the catalog. The boxed 
and labeled areas enclose coherent structures in the data cube. Objects outside these areas also belong most likely to the LA. Note: 
White colors inside regions with disk emission denote column densities in excess of Nui = 9 • 10 22 cm -2 . 



Performing a mass-weighted, least square approximation of 
the vgsr(^ms), we find the relation accounting for all LA objects: 

vgsr(/ms) = 84.2 km s" 1 - 1.2 km s" 1 • / M s. (1) 

The linear regression line is shown in Fig. [2] (bottom panel). 
Despite the high number of low-mass objects, the high-mass 
objects dominate the linear approximation of the overall velocity 
distribution of the LA. Considering the global motion of the LA, 
the velocity at the origin of the LA, with vgsr(0) = 84.2 km s" 1 
is consistent with the value vgsr = 84 + 7 km s" 1 as derived by 
van der Marel et al. 1 d2002h from an ensemble of carbon stars of 
the LMC. This overall velocity structure of the Hi gas suggests 
that the LA has its origin in the LMC rather than in the SMC 
(Put man et al.ll 19981) . This is consistent with the morphological 
and positional c oincidences of the L AI feature with the LMC, 
as suggested by iNidever et al.l (120081) . who have used quite dif- 
ferent methods for their analysis. 

Quantitatively comparing the velocity structure of the 
gaseous objec t s of the LA with the N-body simulation of 
IConnors etaD (l2006h (their Fig. 9) - who attribute the LA to 
the SMC - shows that the simulated velocity distribution is in 



excess of about 200 km s" 1 from the observed one. Their LMC 
orbit is, however, consistent with the Hi velocity structure of the 
high-mass objects. 

4. 1 . Velocity structure of LA 1 

We distinguished three LAI subpopulations as shown in 
Fig.p] (right panel). Toward s a single filament of LA 1.3 
McClure-Griffiths et al. ( 2008) propose that this filament inter- 
acts with the Milky Way galaxy disk. Hi mass-weighted linear 
fits the individual LA filaments underline the general trend of 
decreasing vgsr with increasing /ms- As is obvious in Fig.[5](top 
right panel) LA 1.3 shows the strongest velocity gradient of all 
LA filaments. Based on our statistical analysis of the LAI. 3, the 
fila ment as a whole appears to be s trongly decelerated consistent 
with lMcClure-Griffiths et aD (l2008h . 

Considering the overall velocity gradient of the LA as an 
ensemble of objects (Fig. [2]) and the large angular separation 
between the LA features (Fig.[T]), we propose that the LA is 
most likely a huge coherent structure. LA 1.3 appears to pass 
through sufficiently dense gas of the Milky Way galaxy. LA 1.2 
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linear fit o LA1.3#1285 




70° 60° 50° 40° 30° 20° 10° 0° 

Magellanic longitude 

Fig. 2. Top: Position (/ms) versus velocity (vqsr) diagram showing the maximum brightness temperature distribution across 
-42.5° < Z?lmc ^ 42.5° across the whole velocity range of interest (starting at T B - 200 mK up to 3 K on a linear scale). On 
the righthand side the LMC Hi emission, while the Milky Way emission is visible towards the bottom. The narrowly confined high- 
velocity dispersion emission lines are associated with nearby galaxies covered by GASS. Obvious is a continuous velocity gradient 
of the Hi emission above the Galactic plane from the LMC systemic velocity to vqsr = km s" 1 . The location of clouds belonging 
to the Wannier complexes WA, WB, WC, and WD are annotated. Bottom: The colored circles denote the location of the individual 
LA catalog objects while the gray- scale figure represents the "edge objects", the Milky Way and the LMC. Each color marks an 
individual LA cloud complex (the ID of the catalog - available at the CDS - entry is given), aiming to demonstrate the coherent 
velocity structure of the complexes. The circle diameter scales according to d oc M~ 3 . The regression line marks the mass-weighted 
velocity gradient with its LMC value of vqsr = 84.2 ± 7 km s" 1 demonstrating that the LA gas most likely is released from the LMC 
gaseous body. 



and LA 1.1 are apparently located at larger z-distances from the 
Galactic plane. The same is true for LA2 and LA3, but now on 
the opposite side of the Galactic disk, as evident from Fig. [3] 

These filaments must have already passed the equator of the 
Galactic coordinate system, but we do not know the galactocen- 
tric distance. In the case o f an exponential radial d ensity distri- 
bution of the Galactic disk ( Kalbe rla" & Kerpll2009l) , the strength 
of the interaction depends critically on the radial distance from 
the Galactic center. Our findings (see Fig.[2]and Sect.© are con- 
sistent with the hypothesis that LA2 and LA3 filaments proba- 
bly passed through a thinner volume density portion at longer 
galactocentric distances than LA 1.3. Using the Kalberla & Kerp 
( 2009) model, we suggest for LA2 and LA3 a galactocentric dis- 
tance of > 40kpc where a density n < 10" 4 cm -3 is expected. 

4.2. Age and size of the whole LA structure 

Figure[T] (right panel) shows a surprisingly high degree of sym- 
metry of the Hi object distribution. LA4 covers the same angular 
range in |Z?msI but is located at the opposite side of the LMC 
as LAI. To scatter the Hi gas across more than 60° on the sky, 
symmetrically distributed around the LMC orbit (Z?ms = 0°), 
one needs a fast acceleration mechanism to overcome the grav- 
itational attraction of the MCS and to widely distribute the 



gas across tens of kpc . Assuming a total gravitational mass of 
M M cs = 2 • 10 10 M o (Ivan der Marel et al.|[2002b the escape ve- 
locity is VLMC,esc - 280 km s" 1 . 

The high degree of symmetry in radial velocity between LA 1 
and LA 4 implies that the whole LA is viewed nearly face on. 
Evaluating the ratio of the radial to tangential velocity vec- 
tor compo nent with the escape veloc ity |vLMC,tanl + |vLMC,radl ^ 

I^LMC,esc 

I (Ivan der Marel et al.l (l2002h gives v L MC,rad = 84 ± 
7kms _1 and VLMC,tan = 281 ± 41kms _1 ), and we can sup- 
port this visual impression quantitatively. Almost all of the to- 
tal velocity vector is in the tangential component. If we assume 
that we watch a nearly "face-on LA system" with an inclina- 
tion against the sky of only 13.6°, we can estimate the mini- 
mum time needed for an LA cloud to travel from the LMC body 
to its current location. For this, we assume a) that the linear 
separation between LA 4 and LAI is less than 52 kpc, and b) 
we neglect drag forces and assume a constant objects orbit ve- 
locity equal to the escape velocity. This yields t = 6 • 10 7 a. 
According to this picture LA 1.3 has a radial separation from the 
Galactic center of 19.3 kpc, while LA4 lies at a galactocentric 
distance of 71. 9 kpc. The LAI. 3 distance estim ate we derived is 
consistent with lMcClure- Griffiths et al.l (120081) . who determined 
R ~ 17 kpc via a kinematic distance estimate. 
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Galactic latitude 




Magellanic longitude 

Fig. 3. This color-coded figure illustrates the applied method to differentiate between Milky Way, LMC, and LA Hi emission. 
The cataloged objects show up with significant different GSR velocity distributions (color), allowing them to be separated from the 
neighboring Hi emission within the data cube. Applying the strong criterion to separate between Milky Way and LA emission (FigQ] 
right) we might miss coherent structures, i.e. /ms = 35° and Z?ms = 28°. These objects are continuously connected in Hi emission 
with the Milky Way galaxy. 



4.3. Physical conditions of the LA4 clouds 

4.3.1 . Evaluating the LA4 ambient volume density 

An estimated galactocentric distance of 72 kpc for LA4 corre- 
sponds to a radial distance of 74 kpc and an altitude tan(20°) x 
74 kpc ^ 26.5 kpc above the disk. Figure[5](bottom left panel) 
suggests that even the faint and distant LA 4 clouds are "deceler- 
ated" by drag forces. This implies that even at such large distance 
above the disk a sufficient high volume density exists to cause 
drag. This hypothesis is supported by individual HT clouds as- 
sociated with LA4. Figured shows the Hi 21 -cm column density 
distribution of two clouds of LA4, where apparently both appear 



to be associated with a cometary shape. They appear to "travel" 
in parallel and consist of fast-moving heads and trailing tails, 
implying that gas is stripped off the clouds, heated up, and con- 
centrated within the cloud's slipstream. This gives us the chance 
to estimate the physical conditions at those long distances from 
the galaxy, so we derive from the angular extent of the head of 
0.5°, with 74 kpc distance, A/hi = 5.2 • 10 19 cm -2 , a volume den- 
sity of ?2hi = 0.03 cm" 3 . The mass of the HT clouds is about 
M - 4.5 • 10 5 M . We estimate the kinetic temperature from the 

line width (FWHM = 6.6 km s" 1 ) to 7^ < * 1000 K. 

This yields a gas pressure of P = 60Kcm~ 3 in the certainly 
partly ionized gas. To constrain the location of the HT features 
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Velocity: +84.10 km/s 




LON — CAR 

Fig. 4. Hi brightness temperature distribution of two HT clouds associated with the LA4 are displayed in gray-scale on a linear 
scale. Both clouds show a cometary shape in this velocity channel of 84 km s" 1 . The contour lines mark Hi gas line temperatures at 
100, 200, and 400 mK. The closest bow-like object connected with the head of one HT cloud is at vqsr = 144 km s" 1 , and all the 
more distant structures are at vgsr = 25 km s _1 .The solid line marks the orbital motion of the MCS, the long dashed line the lagging 
rotating halo, and the short dashed line the projection of the resulting velocity vector. 



we assume that they move with the same velocity as the MCS. 
The most probable position along the line of sight would be that 
closest to the MCS orbit. We obtain a distance of ~ 4 kpc with 
a halo density of nuaio = 3 • 10" 3 cm -3 dKalbe rla 2003). This den- 
sity would be just high enough to cause significant drag forces 
(lOuilis & Moord l2Q0lh . FigureH shows that the Z MS = 23°, 
^ms = 30° cloud is positionally correlated with a bow shaped 
feature at vgsr = 144 km s" 1 (contour lines), suggestive of an 
interaction. Accordingly, the objects might interact with local 
density enhancements. An even larger bow-like emission, fea- 
tures are located ahead at velocities vqsr = 25 km s" 1 . 

The orientation of the projected HT velocity vector indi- 
cates its motion with respect to the ambient medium.The result- 
ing vector is aligned with the tail of the HT feature. For this 
estimate we use t he me an 3-D velocity vector determined by 
iKallivavalil etaD ([2006) with |v| = 378 km s _1 . 



5. Head-tail clouds towards the LA 

Head-tail clouds are asymmetric structures. Such asymmetries 
are usually determined by visual inspection, but such a descrip- 
tion may contain personal biases. As an attempt to determine 
the asymmetry of these structures in an objective way, we in- 
cluded higher moments in our analysis. The skewness in partic- 
ular appears to be an appropriate measure for head-tail clouds. 



A well known problem, however, is that higher moments are in- 
creasingly stronger affected by uncertainties. Unfortunately, the 
signal-to-noise ratio is low for most of our clouds. Skewness 
is then ill-defined. The column density of the tail is typically 
lower than that of the head. For a lower signal-to-noise ratio it 
gets increasingly harder to measure asymmetry since the column 
density of the tail may eventually be undetectable. We searched 
for another more robust measure of asymmetry using parameters 
that are less affected by noise. Using our catalog of LA clouds, 
it is feasible to systematically search for clouds that disclose a 
positional difference between the on the sky projected center of 
mass (CM) and the brightness temperature maximum (MI, max. 
intensity). Both parameters are not severely affected by noise but 
a displacement represents a measure for asymmetries. 

In total we count 136 clouds with angular sizes in excess of 
14.4 arcmin and with |CM - MI| > 14.4arcmin. This implies 
that about 30% of our cataloged objects can be considered mor- 
phologically as HT objects. 

This rate is comparable to iPutman et al.l (1201 ll) who have 
found that 35% of the HIPASS compact and semi-compact 
HVCs can be classified as HT clouds from their morphology. 
However, their 116 HT objects were selected from a total num- 
ber of about 2000 HVCs on the whole southern sky. Our sam- 
ple of 136 HT objects was derived from a region covering 1.6 
sr, only about one quarter of the area studied by Put man et al.l 
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Fig. 5. Velocity gradients of the individual coherent structures (see FigO. Top left: LA1.1 and LAI. 2 are located close to the 
Magellanic Clouds. Their ensemble properties are unique, and they dislose a rising vgsr velocity with increasing /ms- Top right: 
LAI. 3 shows the strongest velocity gradient, consistent with the impact hypothesis proposed by iMcClure- Griffiths et al.l (T2008). 
Middle panel: LA2 (left panel) and LA3 (right panel) are localized at large Magellanic stream longitudes. Bottom left: The newly 
classified feature LA4 shows a velocity gradient comparable to LA2 and LA3 despite its close location to the Magellanic Clouds. 
Bottom right: Those cataloged objects not accounted to the coherent structures defined in Fig. Q] disclose the same overall velocity 
gradient implying their association with the LA. 



(2011). We conclude that the LA region has a significant excess 
of HT structures. 

The fainter one of the two HT features shown in Fig. [4] has 
remained undetected so far. 



H T HVCs are comm on structures of HVCs (iBriins et al.l 

2000; Putman et al.ll201ll) . Their cometary shape implies a phys- 
ical interaction of the HVC gas wit h the dynamical different 
ambient medium. Bri ms et al T (120001) show that a linear rela- 
tion exists between the peak column density of HVCs and the 
fraction of observed HT HVCs. This implies that the higher the 
peak column density, the higher the probability of identifing an 
HT structure. This finding suggests that the HVC gas is com- 
pressed at the head, the gas cools more efficiently due to the 
higher volume density, and this leads to lower kinetic temper- 
atures (narrower line widths) and apparently highe r Hi bright- 
ness temperatures. Only towards HVC 125+41-207 (iBriins et al.l 
2001) do we pres e ntly h ave evidence of this process in great 
detail. IBriins et al.l (120011) show the proposed temperature struc- 
ture using a Gaussian decomposition of the Effelsberg Hi data. 
On a statistical basis, the cometary shape on its own is no more 
than a morphological implication, nor does the linear correlation 
between peak column density and HT fraction directly disclose 
ram pressure as a physical process. 



To search for physical indicators for ram pressure interaction 
with the ambient medium, we investigated the ratio between the 
brightness temperature (7b) and the kinetic temperature (T^ n ). 



T B = r s -(l-0<7s = 7ki n 
mn • Av 2 
~ 8kln(2) 



(2) 



Here T$ denotes the spin temperature, the Boltzmann popula- 
tion of the two energy levels, Av is the FWHM of the Hi emis- 
sion line, and r the opacity of the Hi gas. Tb can be considered 
accordingly as lower, while T^ n is an upper temperature limit 
of the HT cloud. The ratio of both temperatures is used as a 
marker to differentiate between high volume density portions of 
the clouds (T^ dominant) and warmer low volume density por- 
tions (Tkin dominant). Using the LA catalog, we calculated maps 
of the temperature ratio across the whole LA region of inter- 
est. We focused in particular on LA 1. 3, because of the com- 
pelli ng morphological argumentation of IMcClure- Griffiths et al.l 
(2008). Figure[6] shows the temperature ratio. Obviously the HT 
cloud (/ms = 27°, Z?ms = -22°) shows a suggestive cometary 
shape and a "cool" head associated with a trailing "warm" tail. 
"Warm" rims are obviously common in LA 1.3 and consistently 
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Fig. 6. Ratio of as a measure of temperature and density gradients. For comparison of the linear scale we plotted the Parkes 
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degree of small-scale structure. The high column density objects are associated with cool cores (heads) localized in the direction of 
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oriented towards the Milky Way disk, indicating physically the 
interaction of the LA gas with the Milky Way disk matter. 

Extending this investigation to the whole field of interest 
(Fig-El discloses a surprising richness of temperature structure. 
All HT structures have warm rims. Objects with higher masses 
have a cold head embedded within warm gas, followed by an 
even warmer trailing gas component. The low-mass objects do 
not show cold condensations. 

Finally we studied the sky parallel component of the vector 
connecting CM and MI and found that the high mass objects 
disclose a coherent velocity pattern (Fig. [7]). The HT velocity 
vectors trace the orbital motion of the LA as a whole. In com- 
bination with the derived overall velocity gradient, the velocity 
segregation of the object masses and the warm rims of the HT 
objects, we conclude that the whole LA is decelerated by drag 
forces towards the Milky Way disk. Both LA 2 and LA 3 have al- 
ready passed the Milky Way disk (see Fig.0. LA 4 fits perfectly, 
spatially as well as kinematically, in the whole picture and ex- 



tends the LA to a huge coherent structure across the southern 
hemisphere. 

6. Summary and conclusions 

We used the GASS data (iMcClure-Griffiths et al.1 120091: 
Kalberla etaD|2010|) to investigate the LA structure. The ana- 
lyzed area covers 70° x 85° degrees. Using a sophisticated image 
characterization tool, it was feasible to differentiate between the 
diffuse Milky Way emission and the LA structures. In total we 
identified 449 objects above the 4 - cr threshold of the original 
data (7b = 200 mK). The applied method not only allows to de- 
termine the location of the objects, but also compiles derived in- 
formation such as radial velocity, phase-space volume, column 
density, geometric center, and temperature weighted center. This 
allowed us to perform an automated statistical analysis of the 
LA as a whole and individual substructures down to the scales 
of individual clouds (HT structures). 
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Fig. 7. Ratio of 7b and T^m across the whole LA as shown in detail in Fig. [6] The arrows mark the orientation of the HT structures 
from the center of mass to the maximum brightness temperature. 



Next to the well studied LAI, LA2, and LA3, we discov- 
ered LA4, which is localized at high Magellanic latitudes but at 
comparable longitudes to LAI. This symmetry but also the asso- 
ciated velocity gradient, suggests that this feature belongs to the 
LA complex. LA4 on average contains less extended clouds. In 
Sect. [4] we argued for a distance between 40 and 70 kpc, which 
implies that these clouds could also have a lower average mass. 

Averaging across the whole LAI to LA4 cataloged clouds, 
we derived a clear negative velocity gradient (deceleration) with 
Magellanic stream longitude. Assuming D = 50 kpc, the mass 
weighted linear correlation yields vqsr = 84.2 km s" 1 as axis 
intercept. This value is identical with in its uncertainties with the 
bulk velocity of the LMC derived by Ivan der Marel et al. I d2002h 
from the carbon star ensemble. According to this we support the 



hypothesis of iNidever et al.l (120081) that the LA as a whole has 
its origin in the LMC. 

The mass of the LA is abou t a fac tor of 50% higher than 
previously a ssume d dBrtins et al.l (120051) at a distance of 50 kpc, 
ICioni et all (120001) ). The main reason for this increase is that 
our survey covered a significantly larger region than is avail- 
able for previous, more specialized investigations. Some objects, 
blended with Milky Way emission features, may be missing in 
our analysis. The LA mass of 38 • 10 6 M derived here is proba- 
bly a lower limit only. 

The velocity structure of individual LA filaments can be in- 
terpreted as prototypical for deceleration by drag forces. LA 1.3 
shows t he strongest velocity gradient, consistent with the pro- 
posal of iMcClure-Griffiths et al.l (120081) . All coherent structures 
of the LA but LA1.1 and LAI. 2 show a systematic deceleration 
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towards the Galactic plane. We find an obvious difference in the 
mass vs. velocity distribution for LAI and LA4 in comparison 
to LA2 and LA3. The last two show that the higher the mass, 
the higher the vqsr velocity. In the framework of drag models, 
these high-mass compact objects keep their momentum, while 
the low-mass, more diffuse clouds are significantly decelerated. 

The apparent symmetry of the LA structure is a surprising 
finding. The whole structure covers nearly 65° across the whole 
sky. At a distance of 50kpc, this corresponds to a linear scale 
of 52.5 kpc. Combining the radi al velocity information o f the 
Hi data with the proper motion (Ivan der Marel et af]|2QQ2h . we 
estimate an inclination of only 13.6° against the plane of the 
sky. With the escape velocity for the LMC we can calculate that 
the LA structures are at least 6 • 10 7 a old. This nearly face-on 
geometry of the LA localized LA 1.3 at a galactocentric distance 
of R ~ 19.3 kpc and LA4 at 7? ~ 71.9 kpc. 

We find pronounced HT structures associated with LA4. 
Using a distance of 74 kpc, an estimate on the kinetic tem- 
perature by the Hi line width, we estima te a pressure of P = 
60Kcm . According to Kalberla (2003) the ambient volume 
density with n ^ 10" 3 cm -3 is high enough to cause drag. 

The analysis of the ratio between and T^m implies that 
the gas of the objects is compressed towards the leading edge 
of the objects. Warm rims are at the head of the HT structures, 
suggesting gas compression and eventually an enhanced cooling. 
All massive LA objects show HT structures. The sky projected 
vector component between the location of their centers of mass 
relative to the maximum Hi line intensities are all oriented par- 
allel to the orbital motion of the MCS. In combination with the 
velocity segregation, the overall velocity gradient and the cool 
leading rims of the HT objects, we conclude that the whole LA 
is a coherent structure decelerated by drag. The newly identi- 
fied LA 4 fits in the whole structure and physical interpretation 
perfectly. 

The emerging picture of a four-arm structure, with two sym- 
metric arms approaching the Milky Way disk (LAI and LA4) 
and two arms that have already passed the disk (LA2 and LA3) 
provides new constrain ts for simulations of t he orbital history of 
the Magellanic Clouds. [Diaz & Bekkil (120111) present a new tidal 
model in which, for the first time, structures resemble a clear bi- 
furcation of the LA. Concerning LA4, we wo uld like to note tha t 
the existence of this feature was predicted by lBesla et al.1 (l2010h . 
but their simulations fail to recover LAI. Hydrodynamic pro- 
cesses such as ram pressure are needed to shape the Magellanic 
Stream. 
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Table 2. Excerpt of the LA catalog of objects available at the CDS. The table comprises several parameters that describe the 449 LA objects 
discussed. 



Units 



Label 



Explanations 



Object 1 



Object 2 



deg 2 km s" 1 

solMass 

sr 

pc z 

deg 

P °_ 
cm -2 

cm -2 

cm -3 

Kern" 3 



Seq 
ID 

Name 

PhaVol 

MassMsol 

SolAng 

Area 

SphDiam 

SphDiampc 

PeakNHI 

MeanNHI 

Density 

Press 

LocMax 



Sequential number 

Internal ID number 

High- velocity cloud (HVC) name (1) 

Phase space volume in deg 2 km s" 1 

Mass at 50 kpc in solar masses 

Solid angle 

Area at 50 kpc from solid angle 
Spherical diameter from solid angle 
Spherical diameter at 50 kpc from area 
Peak column density 
Mean column density 
Mean cloud density at 50 kpc 
Mean cloud pressure at 50 kpc in Kcm -3 
Local maxima > 0.24 deg and > 2.4 km s" 1 



1 

258 

HVC 291.7+21.5+111.1 
0.11 
323.39 
2.07E-05 
51795.3 
0.29 
256.8 
1.53E+018 
7.79E+017 
9.80E-04 
0.05 
3 



2 

259 

HVC 293.4+37.1+086.4 
1.19 
3848.43 
5.23E-05 
130910.9 
0.46 
408.2 
7.81E+018 
3.66E+018 
3.00E-03 
5.23 
6 



- 


RegMax 


Regional maxima > 0.56 deg and > 5.6 km s" 1 


1 


2 


deg 


MSLON 


CM Magellanic stream longitude 


53.04 


69.30 


deg 


MSLAT 


CM Magellanic stream latitude 


-14.93 


-16.40 


deg 


GLON 


CM Galactic longitude 


291.65 


293.42 


deg 


GLAT 


CM Galactic latitude 


21.48 


37.13 


deg 


RAdeg 


CM right ascension (J2000) 


179.08 


184.49 


deg 


DEdeg 


CM declination (J2000) 


-40.17 


-25.12 


kms" 1 


HRV 


CM heliocentric radial velocity 


115.81 


88.05 


kms -1 


LRV 


CM local standard of rest velocity 


111.08 


86.35 


kms -1 


GRV 


CM Galactic standard of rest velocity 


-79.17 


-74.58 


kms -1 


LGRV 


CM local group standard of rest velocity 


-147.88 


-144.62 


kms -1 


FWHM 


CM velocity FWHM 


1.64 


11.54 


- 


xxVar 


CM xx Variance (2) 


0.17 


1.03 




yyVar 


CM yy Variance (2) 


1.16 


1.20 




xyVar 


CM xy Variance (2) 


0.19 


-0.19 




xSkew 


CM xSkewness (2) 


0.57 


-0.00 




ySkew 


CM ySkewness (2) 


-0.23 


-0.08 




xKurt 


CM xKurtosis (2) 


-0.97 


-0.87 




yKurt 


CM yKurtosis (2) 


-1.09 


-1.01 


deg 


PA 


CM position angle of ellipse to north (3) 


150.6 


38.9 




Ecc 


CM eccentricity of ellipse (3) 


0.94 


0.62 


deg 


MinDiam 


CM minor diameter of ellipse (3) 


0.16 


0.41 


deg 


MajDiam 


CM major diameter of ellipse (3) 


0.46 


0.53 


deg 


PTBLON 


Peak T B Magellanic stream longitude 


53.00 


69.24 


deg 


PTBLAT 


Peak Tb Magellanic stream latitude 


-14.96 


-16.40 


kms -1 


PTBGRV 


Peak T B GSR velocity 


-78.32 


-73.38 


deg 


PTBAng 


Angle between CM - peak Tb line and north 


119.5 


90.0 


K 


PeakTB 


Peak brightness temperature (T B ) 


0.32 


0.45 


kms" 1 


PTBFWHM 


Peak T B velocity FWHM 


1.64 


9.07 


K 


TDopp 


Peak Tb Doppler temperature at 50 kpc 


59.3 


1796.3 


pc 


SizeLWRel 


Peak Tb Size-linewidth relation 


1.6 


287.7 


deg 


CenLON 


Centroid Magellanic stream longitude (4) 


53.04 


69.30 


deg 


CenLAT 


Centroid Magellanic stream latitude (4) 


-14.93 


-16.40 


kms -1 


CenGRV 


Centroid GSR velocity (4) 


-79.18 


-74.57 


deg 


SizeLON 


Size along Magellanic stream longitude 


0.23 


0.61 


deg 


SizeLAT 


Size along Magellanic stream latitude 


0.48 


0.40 


kms -1 


SizeVel 


Size along radial velocity 


3.29 


12.36 
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